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Dobrava  virus,  recently  isolated  from  a  yellow- 
neck  mouse  (Apodemus  flavicollis),  captured  in 
a  northern  Slovenian  village  where  severe  cases 
of  hemorrhagic  fever  with  renal  syndrome  were 
recognized,  was  shown  by  serology  and  restric¬ 
tion  enzyme  digestion  of  PCR-amplified  gene 
segments  to  be  related  to  previously  recognized 
hantaviruses.  To  investigate  further  the  relation¬ 
ship  of  this  new  isolate  to  other  hantaviruses,  a 
portion  of  the  medium  (M)  genome  segment  of 
Dobrava  virus  was  amplified  by  PCR  and  the  nu¬ 
cleotide  sequence  determined.  Comparing  the 
nucleotide  sequence  with  the  same  gene  region 
of  other  hantaviruses  revealed  an  overall  homol¬ 
ogy  of  41.7%.  A  phylogenetic  tree  based  on  pair¬ 
wise  sequence  similarity  clearly  showed  that  Do¬ 
brava  virus  is  genetically  distinct,  and  probably 
represents  a  new  virus  in  the  genus  Hantavirus 
of  the  family  Bunyaviridae.  ©  1993  Wiley-iiss,  Inc 
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INTRODUCTION 

The  genus  Hantavirus  is  comprised  of  at  least  6  ge¬ 
netically  and  antigenically  related  but  distinct  viruses: 
Hantaan  (HTN),  Seoul  (SEO),  Puumala  (PUU),  Pros¬ 
pect  Hill  (PH),  Thailand  (THAI),  and  Thottapalayam 
^)(THP)  [Lee  et  al.,  1985;  Schmaljohn  et  al.,  1985;  Xiao 
et  al.,  1992).  Some  of  these  viruses  are  etiologic  agents 
of  the  human  disease,  hemorrhagic  fever  with  renal 
.^syndrome  (HFRS).  Recently,  a  new  isolate,  designated 
•'^Dobrava  virus,  was  obtained  from  an  yellow-necked 
mouse,  Apodemus  flavicollis,  captured  in  Slovenia 
_^[Avsic-Zupanc  et  al.,  1992],  where  severe  cases  of  HFRS 
pjj were  identified.  Preliminary  serological  studies  and 
^polymerase  chain  reaction-restriction  fragment  length 
-^^polymorphism  (PCR-RFLP)  analysis  demonstrated  its 
relationship  with  the  existing  hantaviruses,  suggest¬ 
ing  that  it  might  represent  a  new  member  of  the  genus 
Hantavirus,  family  Bunyaviridae  [Avsic-Zupanc  et  al., 
1992). 


We  now  report  the  nucleotide  sequence  comparison  of 
a  PCR-amplified  fragment  from  the  M  genomic  seg¬ 
ment  of  Dobrava  virus  with  previously  published  se¬ 
quences  of  other  hantaviruses.  Dobrava  virus  has  found 
to  be  genetically  distinct  from  other  hantaviruses,  add¬ 
ing  further  evidence  that  Dobrava  represents  a  new 
hantavirus. 

MATERIALS  AND  METHODS 
Virus  and  Cell  Culture 

Dobrava  virus  was  isolated  from  the  lungs  of  an  Apo¬ 
demus  flavicollis  [Avsic-Zupanc  et  al.,  1992).  The  origi¬ 
nal  specimens  were  collected  from  Dobrava  village  in 
Slovenia,  where  severe  cases  of  HFRS  were  identified. 
Vero  E6  cells  were  infected  with  the  virus  and  pro¬ 
cessed  for  total  cellular  RNA  extraction  by  using  a 
guanidine  isothiocyanate/cesium  chloride  method,  as 
described  previously  [Xiao  et  al.,  1991). 

Reverse  Transcriptase  Polymerase  Chain 
Reaction  (RT-PCR)  and  Direct  Sequencing 
of  the  Amplimer 

Two  primer  sets  used  in  the  initial  PCR  amplifica¬ 
tion,  HanG2-l:2  and  HanG2-3:4,  were  derived  from  the 
glycoprotein  G2  encoding  region  of  the  Hantaan  virus 
M  genome  segment  [Xiao  et  al.,  1991).  The  sizes  of  the 
expected  amplimers  were  411  and  334  bp,  respectively. 
RT-PCR  reactions  were  done  as  described  previously 
[Xiao  et  al.,  1991).  Briefly,  total  cellular  RNA  was  cop¬ 
ied  to  cDNA  with  avian  myeloblastosis  virus  reverse 
transcriptase,  then  amplified  for  35  cycles  with  the  Taq 
polymerase,  in  an  automatic  DNA  Thermal  Cycler 
(Perkin-Elmer-Cetus,  Norwalk,  CT). 

The  amplified  product  was  subjected  to  preparative 
electrophoresis  by  using  1 .2%  low  melting  point  (LMP) 
agarose  gel  (IBI,  New  Haven,  CT)  prepared  with  Tris- 
acetate-EDTA  (TAE)  buffer.  The  target  band  was  cut 
out  and  DNA  recovered  by  using  the  GENECLEAN  Kit 
(BIO  101,  La  Jolla,  CA).  Ten  microliters  of  double- 
stranded  DNA  (1.5  u-g)  was  denatured  with  2  p.1  of  2  N 
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Fig.  1.  Identification  of  cognate  region  for  sequence  comparison  and  dendrogram  construction  Se¬ 
quences  were  aligned  against  that  of  Hantaan  7tf  - 1 1  h  <  Hash  value  6i 


NaOH,  ethanol  precipitated,  and  dried  briefly.  Se¬ 
quencing  was  done  by  the  dideoxy  chain  termination 
method  1  Sanger  et  al.,  19771,  with  the  Sequenase  Ver¬ 
sion  2.0  DNA  Sequencing  Kit  from  United  States  Bio¬ 
chemical  (Cleveland,  OH).  Primers  HanG2-l  and 
HanG2-2  were  used  for  the  sequencing  reactions. 

Sequence  Analysis 

Reference  sequences  were  previously  published  and 
deposited  in  GenBank.  Sequences  examined  included: 
Hantaan  virus  strains  76-118,  Lee,  and  Hojo  ISchmal- 
john  et  al.,  1987,  1988;  Seoul  virus  strains  HR8039  and 
SR-11  lAntic  et  al.,  1991;  Arikawa  et  al.,  1990|;  Puu- 
mala  virus  strains  Sotkamo  and  Hallnas-Bl  IGiebel 
et  al.,  1989;  Vapalahti  et  al.,  19921;  Prospect  Hill  virus 
strain  PHV-I  IParrington  et  al.,  1991 1:  and  Thailand 
virus  strain  Thai749  (Y-K.  Chu  et  al.,  unpublished 
data).  The  Dobrava  virus  sequence  was  compared  with 
all  of  these  by  the  MacVector  computer  program  from 
International  Biotechnologies,  Inc.  (New  Heaven,  CT). 
to  identify  homologous  regions  in  each  sequence.  The 
resulting  sequence  fragments  were  then  compared  with 
each  other  by  using  the  CLUSTAL  program  of  the  PC/ 
Gene  software  from  IntelliGenetics,  Inc.  (Mountain 
View,  CA>.  Pairwise  similarity  scores  were  calculated 
for  all  pairs  of  sequences  and  the  degree  of  similarity 
used  to  form  a  tree  (dendrogram)  by  clustering  the  se¬ 
quences  based  on  overall  similarity. 

RESULTS 

Dobrava  virus  was  successfully  amplified  with  the 
HanG2-l:2,  but  not  with  the  HanG2-3:4  primer  pair, 
while  both  primer  pairs  amplified  Hantaan  virus  strain 
76-118  (data  not  shown).  The  amplimer  generated  with 
HanG2-l:2  was  processed  for  sequencing.  A  sequence  of 
207  bp  from  the  HanG2-l:2  amplimer  was  obtained  and 
aligned  against  the  whole  M  segment  sequence  of  Han¬ 


taan  virus  strain  76-118.  along  with  those  of  8  other 
hantaviruses.  The  cognate  region  was  identified  for 
each  sequence  as  shown  in  Figure  1.  corresponding  to 
the  region  2120  to  2330  in  the  Hantaan  virus  cDNA 
sequence.  This  region  was  extracted  from  all  the  se¬ 
quences  for  further  pairwise  sequence  comparison  by 
the  CLUSTAL  program. 

Alignment  of  the  10  nucleotide  sequences  is  shown  in 
Figure  2.  A  dendrogram  was  constructed  based  on  the 
calculated  scores  of  pairwise  similarity  (Fig.  3'.  The 
length  of  compared  region  was  21 1  bp.  with  an  overall 
identity  of  88  bp  (41 ,7N  >.  To  ensure  the  reproducibility 
of  this  dendrogram  without  Dobrava  virus,  another  500 
bp  region  from  the  glycoprotein  G1  encoding  gene  of  the 
completely  sequenced  strains  was  analyzed  with  the 
same  strategy.  An  identical  dendrogram,  excluding  Do¬ 
brava  virus,  resulted. 

DISCUSSION 

Seroepidemiologica!  investigations  have  identified  at 
least  three  distinct  antibody  reaction  patterns  among 
HFRS  patients  seen  in  Yugoslavia  and  Slovenia  I  Avsic- 
Zupanc  et  al.,  1990;  Gligic  et  al.,  1 992  j.  These  include 
reactions  specific  for  a  virus  similar  or  identical  to  pro¬ 
totype  Hantaan  virus,  a  second  pattern  specific  for  Puu- 
mala  virus,  and  a  third  pattern  that  appears  to  be  Do¬ 
brava  virus-specific.  In  addition,  examining  Dobrava 
virus  with  monoclonal  antibodies  to  various  hantavi¬ 
ruses  demonstrated  a  unique  serological  reactivity  pat¬ 
tern  for  this  virus.  Finally,  PCR-RFLP  comparisons  of 
Dobrava  and  other  recognized  hantaviruses  likewise 
found  that  this  virus  is  related  to,  yet  distinct  from,  all 
previously  recognized  hantaviruses  |  Avsic-Zupanc 
et  al.,  19921. 

The  results  show  that  at  the  genomic  level.  Dobrava 
virus  shared  howmlngv  with  other  hantaviruses  How  - 
ever,  it  was  clearly  distinct  when  the  sequence  of  the 


:S 


A- 


a.o 


DTIC  QUALITY  INSPECTS®  % 


CDE4_-., 


154 


Xiao  et  al. 


DOBRA  VA  CCCATTAATGAAGAACAATCAGTACATCTTCACATAGAGATAAGGAGTCA  50 
HTN  CCACTTGAGGAAGCACAATCCATTGACCTACATATTGAAATAGAAGAACA  5C 
LEE  CCACTTGAAGAAGCACAATCCATTGATCTACATATTGAAATAGAAGAACA  50 

HOJO  CCACTTGAAG - CAC AATCC  ATTGATCT ACAT  ATTG AAAT  AGAAGAACA  47 

SRI  1  CCAGTTAATGACCAACAGAGTGTCTCATTGCATATAGAAATTGAAAGTCA  50 
SEOUL  CCAGTTAATGACCAACAGAGTGTCTCATTGCACATAGAAATTGAAAGTCA  50 
THAI749  CCAGTTAACACAGAACAAGGAGTTCAATTACATATAGAGATTGAAGAGCA  50 
HALLNAS  CCAGCAAACGAACAAGAGAAAATCCCATTTCATCTGCAGTTAAGCAAACA  50 
SOTKAMO  CCGGCAAATGAACAAGAAAAAATACCATTTCATTTACAGATAAGTAAACA  50 
PHV  - GCAAATGAGCAAGAAAGAATCCCTTTTCACTTCCAGCTAGAGCGTCA  47 


DOBRAVA  AGGATTTC-TACAAGTGTTCATGCATTAGGTCATTGG7TTGATGGA-GAC  98 
HTN  GACAATTGGTGTTGATGTGCATGCTCTAGGACACTGGTTTGATGGTCGTC  100 
LEE  AACAATTGGTGTAGATGTGCATGCTCTAGGACACTGGTTTGATGGTCGTC  100 
HOJO  AACAATTGGTGTGGATGTACATGCTTTAGGACACTGGTTTGATGGTCGTC  97 
SRI  1  AGGCATTGGTGCTGCTGTTCATCATCTTGGACATTGGTATGATGCAAGAT  1 00 
SEOUL  AGGCATTGGTGCTGATGTCCATCATCTTGGACATTGGTATGATGCAAGAT  100 
THAI7  19  AGGAATTGGTGCAGATGTCCATCACCTAGGCCATTGGTTTGATGCCCGGC  100 
HALLNAS  AGTGATTCATGCAGAGATCCAGCATTTAGGTCATTGGATGGATGCTACAT  100 
SOTKAMO  GGTGATTCATGCTGAAATTCAACATCTAGGCCACTGGATGGATGCTACAT  100 
PHV  GGTGATTCATGCAGAGATTCAAAATTTAGGCCACTGGATGGATGGAACAT  97 


* 


4 


DOBRAVA  TTAATTTAAAGACATCTTTTCATTGTTATGGTGCATGTACTAAGTATGAA  1 4  8 
HTN  TTAACCTTAAAACATCCTTTCACTGTTATGGTGCTTGTACAAAGTATGAA  150 
LEE  TTAACCTTAAAACATCCTTTCACTGTTATGGTGCTTGTACAAACTATGAA  150 
HOJO  TTAACCTTAAAACATCCTTTCACTGTTATGGTGCTTGTACAAAGTATGAA  147 
SRI  1  TGAATCTAAAAACCTCATTTCATTGTTATGGTGCCTGCACAAAATATCAA  150 
SEOUL  TGAATTTAAAAACCTCATTTCATTGTTATGGTGCCTGCACAAAATATCAA  150 
THAI749  TTAATTTAAAAACATCATTCCATTGTTACGGGGCATGCTCAAAGTACCAG  150 
HALLNAS  TTAATCTTAAAACTGCATTTCACTGCTATGGCTCATGTGAGAAGTATGCT  150 
SOTKAMO  TTAACTTAAAAAGTGCTTTTCACTGTTATGGTTCATGTGAAAAATATGCC  150 
PHV  TTAATCTGAAAACATCATTCCATTGTTACGGTGCTTGTGAAAAGTATGCT  147 


DOBRAVA  TATCCTTGGCATACAGCAAAATGCCACTTTGAAAGGGATTTCGAGTATGA  198 
HTN  TACCCTTGGCATACTGCAAAGTGCCATTATGAAAGAGATTACCAATATGA  200 
LEE  TACCCTTGGCATACTGCAAAGTGCCACTATGAAAGGGATTACCAATATGA  200 
HOJO  TACCCTTGGCATACTGCAAAGTGCCACTATGAAAGGGATTACCAATATGA  1 97 
SRI  1  TACCCATGGCACACTGCAAAATGCCATTTTGAGAAAGATTATGAGTATGA  20  0 
SEOUL  TATCCATGGCACACTGCAAAATGCCATTTTGAGAAAGATTATGAGTATGA  200 
THAX749  TACCCTTGGCATACAGCCAAGTGTCACTTTGAAAAAGATTATGAATATGA  200 
HALLNAS  TATCCTTGGCAGACAGCAGGTTGTTTCATAGAAAAAGATTATGAATATGA  200 
SOTKAMO  TACCCCTGGCAAACAGCAGGATGTTTCGTTGAAAAAGATTATGAATATGA  200 
PHV  TATCCGTGGCAAACTGCAAAATGTTTCTTAGAGAAAGATTATGAATTTGA  197 


DOBRAVA 

GAACA — TGGG 

207 

HTN 

GACGAGCTGGG 

211 

LEE 

GACAAGCTGGG 

211 

HOJO 

GACAAGCTGGG 

208 

SRI  1 

AAATAGCTGGG 

211 

SEOUL 

AAATAGCTGGG 

211 

THAI749 

AAATAGTTGGG 

211 

HALLNAS 

GACTGGTTGGG 

211 

SOTKAMO 

AACTGGTTGGG 

211 

PHV 

AA^AGGGTGGG 

208 

P  ig.  2.  Multiple  sequence  alignment  of  the  homologous  region.  The  dash  -signs  represent  gaps. 


PCR  amplified  segment  is  compared  with  other  hanta¬ 
viruses.  This  finding,  combined  with  previously  pub¬ 
lished  observations  of  serological  uniqueness,  suggests 
that  Dobrava  is  yet  another  distinct  member  of  the 
gem’s. 

Construction  of  phylogenetic  trees  based  on  nucle¬ 
otide  sequence  similarity  has  been  used  to  depict  rela¬ 


tionships  of  closely  related  viruses  as  an  alternative  to 
the  traditional  serotyping  of  viruses.  Rico-Hesse  et  al. 
|1987|,  comparing  nucleotide  sequences  of  !5G  kp, 
showed  that  this  method  of  comparison  was  able  to 
establish  epidemiologically  consistent  relationship 
among  widely  distributed  isolates  of  wild  poliovirus 
type  1  genotypes.  From  the  dendrogram  constructed  in 
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Fig.  3.  Dendrogram  based  on  the  pairwise  comparison  of  the  hanta¬ 
virus  sequences.  The  horizontal  distance  to  the  node  that  connect  2 
strains  is  proportional  to  the  nucleotide  sequence  divergence  between 


this  study,  Dobra va  virus  was  found  to  be  most  closely 
related  to  isolates  of  Hantaan,  Seoul,  and  Thailand  vi¬ 
ruses;  whereas  it  was  more  distantly  related  to  isolates 
of  Puumala  and  Prospect  Hill  viruses  (Fig.  3).  However, 
we  showed  it  to  be  clearly  separated  from  all  of  these 
viruses,  in  agreement  with  the  serological  results 
[Avsic-Zupanc  et  al.,  1992].  Construction  of  a  second 
dendrogram,  based  on  the  sequence  analysis  of  a  500  bp 
segment  from  the  G1  encoding  region  (excluding  Do- 
brava  virus),  showed  an  identical  relationship  among 
these  viruses,  verifying  the  reliability  of  this  partial 
sequencing  analysis.  In  addition,  similar  relationships 
were  recently  reported  by  others  who  compared  gene 
product  sequences  from  the  whole  M  segment  of  the 
hantaviruses  with  complete  sequence  information  [An¬ 
tic  et  al.,  1992),  further  adding  credence  to  our  proposed 
placement  of  Dobrava  virus. 

The  nucleotide  sequence  is  the  most  reliable  marker 
in  identifying  a  virus;  however,  complete  sequencing  of 
all  new  isolates  is  neither  practical  nor  necessary  for 
this  purpose.  Data  from  this  study  showed  that  the 
overall  interrelationship  of  hantaviruses  documented 
from  a  region  of  approximately  21 1  bp  is  very  similar  (if 
not  identical)  to  that  from  the  whole  M  genome  seg¬ 
ment  of  hantavirus  [Antic  et  al.,  1992].  This  provides  a 
very  useful  and  rapid  tool  for  virus  classification  or 
typing,  with  the  easy  availability  of  purified  DNA  tem¬ 
plate  generated  by  PCR. 
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